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ABSTRACT 

We present a quantitative and relatively model-independent way to assess the radial structure of nearby AGN tori. These putative 
tori have been studied with long-baseline infrared (IR) interferometry, but the spatial scales probed are different for different objects. 
They are at various distances and also have different physical sizes which apparently scale with the luminosity of the central engine. 
Here we look at interferometric visibilities as a function of spatial scales normalized by the size of the inner torus radius R m . This 
approximately eliminates luminosity and distance dependence and, thus, provides a way to uniformly view the visibilities observed 
for various objects and at different wavelengths. We can construct a composite visibility curve over a large range of spatial scales if 
different tori share a common radial structure. The currently available observations do suggest model-independently a common radial 
surface brightness distribution in the mid-IR that is roughly of a power-law form r" 1 as a function of radius r, and extends to ~100 
times R m . Taking into account the temperature decrease toward outer radii with a simple torus model, this corresponds to the radial 
surface density distribution of dusty material directly illuminated by the central engine roughly in the range between r° and r . This 
should be tested with further data. 
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1. Introduction 

Speckle and long-baseline interferometry in the infrared (IR) 
has started to explore the spatial stru cture of the putative tori 
in Active Galactic Nuclei (AGN; e . g. IWitlkowski et all 1 1998t. 
IWeinberger et allll999l;ISwain et all l2003t I Weigelt et all 12004; 
iJaffe et all l2004t IWittkowski etaiTl2004t iMeisenheimer et al.L 
2007; rTristram et alll2007l) . However, the number of long base- 
lines per object is still generally limited, and a given baseline 
probes different spatial scales for targets at various distances 
and having different physical sizes. Here we investigate the pos- 
sibility of uniformly view these interferometric measurements 
taken for various objects spanning over different wavelengths, 
by studying them as a function of spatial scales in units of torus 
inner radius. 



2. Composite radial structure 

2.1. Spatial frequency per inner torus radius 

The radius R m of the inner boundary of the torus dust distribu- 
tion is thought to be set by dust sublimation. In this case, for a 
given sublimation temperature and grain size distribution, R m is 
expected to be proportional to L 1 ^ 2 , where L is the U V/optical lu- 
minosity of the central engine (e.g. lBarvainisl 1987b. This propo- 



; propo- 

tionality was recently confirmed by iSuganuma et alJ 12006) 
through near-IR reverberation measurements. These provide the 
inner torus radius as the light travel distance for the time-lag 
between the variability in the optical and the near-IR. It is con- 
ceivable that the dust temperature structure within the torus is 
primarily determined by the illumination from the central source 
and is, thus, also expected to scale with L 1 ' 2 . In this case, if we 
view various interferometric size information as a function of 
spatial scales normalized by R[ n , at least the primary luminosity 



dependence of the data, as well as the distance dependence, is 
eliminated. 

A direct observable in an interferometric measurement is the 
visibility, or the normalized Fourier amplitude of the brightness 
distribution of a source, as a function of spatial scales (or spatial 
wavelength; see below) given in angular size. Here we normalize 
the spatial scale by the angular size of R[ n . More precisely, the 
Fourier component is studied as a function of spatial frequency, 
i.e. the number of spatial cycles in a given angular size. Thus 
we study here visibilities as a function of the number of spatial 
cycles in R\ n , or spatial frequency per R m . 

For a given projected baseline b in m, observing wavelength 
A in fim, and Ri n in milli-arcsecond (mas), the spatial freqency 
in cycles per R m , written here as u ■ R m (where u is a spatial 
frequency per mas), is given as 



u ■ R in = 4.85 x 10" 



i bjm) 
A (fJ.m) 



■ Ri n (mas) (cycles per /?i„). (1) 



The corresponding spatial wavelength A, i.e. the reciprocal of 
spatial frequency u, is given in units of R[ n as 



A//? in = 2.06 x 10 ; 



2 ^ (/™) 



1 



(2) 



b (m) R\ n (mas) 

This corresponds to the representative spatial resolution of the 
configuration (b and A). For simple geometries, the visibility 
takes the first null at a spatial wavelength A roughly equal to the 
characteristic size (e.g. diameter of an uniform disk or a ring). 
The quantity A can also be compared directly to the diffraction 
limit of 2.5 x 10 2 (mas) A(/jm)/D(m) for a single aperture di- 
ameter D. We will study visibilities as a function of these two 
quantities, u ■ R m and A/R[ n , throughout this paper. 

An accurate scale length for the inner boundary radius R[ n 
has not been well determined yet. As the most plausible, obser- 
vationally motivated quantity, we adopt the one obtained from 
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Table 1. List of AGNs used in this paper. 



name type cz" obs. date inst. A b PA R m K-band reference APA S 

(kms~') (UT) (pm) (m) (deg) (mas) Lag (days) (deg) 



NGC1068 


2 


914 


2003-06-15/16 


MIDI 


8.2-13 


79 


-179 


1.6* 




Jaffe et al. (2004) 


89 


NGC1068 


2 


914 


2003-11-09/10 


MIDI 


8.2-13 


46 


-135 


1.6* 




Jaffe et al. (2004) 


45 


Mrkl239 


1 


6321 


2005-12-19 


MIDI 


8.2-13 


41 


36 


0.17 r 




Tristram et al. (2008) 


86 


NGC3783 


1 


3234 


2005-05-28 


MIDI 


8.2-13 


43 


46 


0.32 


85'' 


archive 





NGC3783 


1 


3234 


2005-05-31 


MIDI 


8.2-13 


65 


120 


0.32 


85'' 


Beckert et al. (2008) 


74 


NGC4151 


1 


1242 


2003-05-23 


KI 


2.2 


83 


37 


0.47 


48 f 


Swain et al. (2003) 


50 



"Radial velocity corrected for the cosmic microwave background from NED. 
^Calculated from the estimated intrinsic V-band flux (Pier et al., 1994) and Eq[3] 
^Calcu l ated fr o m the V-band flux f y = 7 8mJv dVeron-Cettv & Veroni 2003) and Eqf3] 
iGlassI <1992l) . iMinezaki etai] j2004l) . ^Suganuma et al.l d2006h . ~ 

s Position angle difference between the projected baseline and the expected torus major axis direction (see text). 



the near-IR reverberation time-lag radii, or, if not available, the 
ov erall fit to them as a fun ction of UV/optical luminosity L given 
by[Suganuma et al. (2006). In the latter case, the angular size for 
the time-lag radius is given as a function of observed flux f v in 
V-band as (the same as Eq.4 in Kishimoto et al. 2007) 

R m = 0.43 [f v (V) I (50 mJy)] 1/2 (mas). (3) 

We note that the time-lag radius is actually a factor of ~3 
smaller, f or a given lumin osity L, than the sublimation radius 
given by iBarvainisI dl987l) assuming graphite grains with ra- 
dial size a=0.05 fim and sublimation temperature r su b=1500 K. 
The possible reasons fo r this difference have been discussed by 
Kish imoto et al.l d2007t) . including the possibility that R m is de- 
termined by much larger grains. As we discuss below, the re- 
verberation radius is supported by the near-IR interferometry of 
NGC 4151 dSwain et alll2003l) . The uncertainty in the time-lag 
radii is ~0.2 dex, based on the scatter of the fit. 

2.2. Data 

The data most suitable for the radial structure study are those 
for Type 1 AGNs. In these objects, the tori are thought to have 
much smaller inclinations than in Type 2 AGNs, so that the ra- 
dial structure can be studied without a significant orientation ef- 
fect. We would also expect that the dependence on the exact PAs 
of the projected baselines is relatively small (the torus image 
projected on to the sky would not be far from being circular- 
symmetric; see more below). Those Type 1 objects for which 
long-baseline IR interferometric data are found in the literature 
or archive are listed in Table [TJ along with the adopted size for 
R m . All the mid-IR data here were obtained with VLTI/MIDI. 
In addition, we also use the M IDI data for the Type 2 AGN 
NGC1068 of lJaffeet all fixM) . The data provide the highest 
spatial frequency information in the mid-IR, and the inclina- 
tion effect at these wavelengths is expected to be relatively small 
compared with that in the near-IR (see more below). 

All the MIDI data we re uniformly re-reduced with the soft- 
ware EWS (version 1.5.2; Jaffe 2004), with a few modifications 
using our own IDL codes. First, since MIDI visibility measure- 
ments for the faint Type 1 targets (< 1 Jy) are largely limited 
by the accuracy of photometric (total flux) measurements, we 
implemented an additional background subtraction for the pho- 
tometry frames. Second, to reduce errors in group delay determi- 
nations, we have smoothed delay tracks over ~10-20 frames. To 
avoid possible positive bias in correlated flux, we also averaged 
~ 10-20 frames for the determination of phase offsets. 

The system visibility was obtained from the observations of 
visibility calibrators usually taken right after or before each tar- 



get observation with a similar airmass. These data were also re- 
duced with similar smoothings above to calibrate out the effect 
of the time averaging. These calibrators ar e also generally th e 
photometric standards found in the list by ICohen et al.l ([1999). 
We obtained the correlated flux, which was corrected for the sys- 
tem visibility, and total flux separately. The errors for the total 
flux and correlated flux were estimated from the fluctuation of 
the measurements over time. For the targets with a few total flux 
measurements available, we took weighted means. The errors 
for the final visibilities were obtained from the estimated errors 
in the correlated flux and total flux (except for NGC 1068 where 
the errors were estimated from the scatter of several available 
visibility measurements). 

2.3. Mid-IR surface brightness distribution 

FigQ] shows the visibility data listed in Table Q] as a function of 
spatial frequency per R^ n (or spatial wavelength A in units of R^; 
upper jc-axis). The x-axis is in log scale to cover a large range of 
spatial frequencies as given by different objects and by different 
observing wavelengths. The red/yellow/green colors correspond 
to wavelenths from 13 to 8.2 yizm, while purple is for 2.2 /mi. The 
MIDI data have been binned with AA ~ 0.4 jjm, with the error 
for each bin taken as the median of the errors over the binned 
spectral channels (since errors over adjacent channels are cor- 
related). For comparison, we plot the visibility curves for ring, 
gaussian, and power-law intensity distributions of different in- 
dices with an inner cut-off radius R m . 

At a given observing wavelength, as the torus inclination in- 
creases from face-on to edge-on cases (i.e. from Type 1 to Type 2 
inclinations), the visibilities are expected to generally decrease, 
since the inner bright radiation becomes obscured and this effec- 
tively makes the overall size of the intensity distribution larger. 
(Here we are restricting our discussion to the spatial frequencies 
within the first lobe, i.e. those less than the one at which the vis- 
ibility reaches the first null.) On the other hand, this decrease 
is expected to be quite small at long wavelengths, such as in 
the mid-IR, due to much smaller obscuration effects. Therefore, 
the mid-IR visibilities of the Type 2 AGN NGC 1068 provide 
quite a meaningful lower limit for those of the Type 1 objects at 
the same spatial frequencies (indicated by the upward arrows in 
Fig©. 

If we put together the lower limit from NGC 1068 with the 
visibilities observed for the Type 1 AGNs in the mid-IR (FigJTJi, 
the data at different spatial freqencies appear quite coherent, 
i.e. consistent with the case where these objects share roughly 
a common radial brightness distribution in face-on inclinations. 
When compared with various model curves shown in FigQ] the 
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Fig. 1. Observed visibilities as a function of spatial frequency 
per inner torus radius. The red/yellow/green colors correspond to 
from 13 to 8.2 pan and purple corresponds to 2.2 fim. The mid-IR 
data are binned with AA ~ 0.4/mi. The arrows are to show that 
the NGC1068 data are taken as the lower limit for the visibilities 
of Type Is at the same spatial frequencies. Visibility curves for 
various simple geometries are also plotted; three dashed curves, 
power-law with index -2.5, -2.0, and -1.5 from top to bottom; 
dotted, Gaussian with HWHM 25 /?;„; dot-dashed, ring with ra- 
dius R m . 



2 object NGC1068 have quite different PAs, which is probably 
adequate for obtaining lower limits for Type 1 visibilities. 

The only existing near-IR measurement for Type 1 AGNs, 
i.e. the one for NGC 4151 (FigQ] purple cross), is consistent with 
a ring of radius 7?; n . It is quite expected that the torus shows a 
much more compact structure in the near-IR than in the mid-IR, 
since the former reflects the distribution of materials at higher 
temperatures. For a given inner boundary radius R\ n , the ring 
model represents essentially the most compact structure, giving 
the upper limit for the visibilities. This essentially means that 
R m generally cannot be much larger than assumed here (the ob- 
served data points would move to the right in Fig Q] if we adopt 
larger for each object, which would exceed the upper limit 
given by the R\ n ring). Future near-IR data at high spatial fre- 
quencies are crucial to settle this issue. 

Note that we should in principle correct the observed visi- 
bility for the possible contribution from the unresolved accre- 
tion disk in the near-IR to estimate the intrinsic visibility of the 
torus alone. The correction depends on the estimate of the near- 
IR flux fraction from the disk, which in turn depends on the as- 
sumed near-IR spectral shape of the accretion disk. The correc- 
tion, ho wever, is estimated to be very small (less than 5% re- 
duction; Kishimoto et al. 2007), as long as the disk spectrum is 
bluer than observed in the UV/optical which seems quite likely 
(e.g. lKishimoto et alj|2008h . 

All the visibility functions shown in FigQ]reach the first null 
at a spatial wavelength A ~ a few R m (outside the figure) due to 
the inner boundary diameter being 2R m . In the mid-IR, detailed 
intensity distributions in the inner few R m region will not affect 
the overall visibility curves significantly, since a major fraction 
of the mid-IR radiation originates from larger regions. To discuss 
whether the same model can account for both the mid-IR and 
near-IR data, we need a physical torus model. 



mid-IR brightness distribution appears to be consistent with a 
power-law form, where the index is ~ -2 at ~8 yum. 

The visibility functions for power-law distributions shown in 
FigQ]are for the cases with the outer cut-off radius R ollt = 100R; n . 
For the power-law distributions with index > -2, R out changes 
the shape of visibility functions, especially at low spatial fre- 
quencies. Here, the data at high spatial frequencies (NGC1068 
and the longer baseline data for NGC3783) roughly fix the power 
law index. Then the low spatial frequency data further constrain 
Rout, which corresponds to the outer radius of the mid-IR emis- 
sion region, to be roughly ~ 1 00i?i n (see also section 3). 

In order to be least sensitive to the possible dependency 
of the visibilities on the position angle (PA) of the projected 
baseline, we would ideally want to compare the visibilities of 
each object as well as different objects measured at the ma- 
jor axis direction of the torus projected on to the sky. The mi- 
nor axis PA can be e stimated from a line ar radio jet structure 
(-7 7° for NGC4151 iMundell et all l200l -0° for NGC1068, 
e.g. iGallimore et alj 120041) or o ptical polarization P A (~136° 
forNGC3783. ISmith et alj |2002h -40° for Mrkl239, [Goodrich] 
119891) . The PA difference between the projected baseline and the 
expected major axis is listed in TableQ] Unfortunately, these rel- 
ative PAs are not uniform for the data gathered here. However, 
they do not seem to disturb the composite visibility diagram 
shown in Fig[T] which might indicate that the PA dependence 
is small as we expect and assumed here for Type Is. Further data 
are needed to address this issue. The two data sets for the Type 



3. Torus model 

We aim here to generically describe the torus surface bright- 
ness distribution with as a few parameters as possible. In 
various torus models, two types of dust distributions have 
been considered; na mely, smooth and clumpy distr ib utions 




(e.g. iPier & Krolikl 1 19921 119931; iGranato & Danese 
lEfstathio~ & Rowan-Rob insonl 119951: iNenkova et alj 
i Dullemond & van Bemmel l2005t [Schartma nn et all 
IHonig et all 120061 ISchartmann et all [2008). We first consider 
the latter and assume that the torus consists of various discrete 
dust clouds each being optically thick. The constraints from the 
alternative case of smooth dust distribution will be considered 
later below. 

The near- and mid-IR torus radiation is primarily dominated 
by that from those c louds which are di rectly illuminated by the 
central engine (e.g. IHonig etail 12006). We can then approxi- 
mately parameterize the face-on surface brightness distribution 
S v (r) at an observing frequency v as a function of radius r as 



S v (r) = I v (T IBax (r))-(r/R io y. 



(4) 



The first term, / y , is the intensity from an illuminated cloud as 
a function of 7' max (r), which is the maximum temperature of the 
dust grains at radius r written as 



(r/RinY 



(5) 



This is derived from the thermal equilibrium of a dust grain and 
involves the wavelength dependency of the absorption efficiency 
of the grain in the form of its spectral index y in the infrared 
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Fig. 2. Observed visibilities as compared with physical torus 
models. Red/green/purple curves are for the observing wave- 
lengh of 13/8.5/ 2.2 //m , respectively. Solid curves are with 
a = 0.0 for the y = case, while dotted curves are with a — -0.9 
for the y — 1.6 case. 



cantly if a larger outer boundary is adopted (the region outside 
100/?i n contains only a small fractional amount in the mid-IR). 
Therefore, in the 100/?j n cases adopted here, the outer bound- 
ary is set almost by the temperature distribution, rather than the 
density distribution of the illuminated material. 

Let us consider briefly the case of smooth dust distribution. 
If the opening angle of the geometrical distribution is constant 
over the radii, i.e. if there is no geometrical flaring, the material 
at large radii is illuminated only by the attenuated radiation from 
the central engine. Then I v would represent the intensity from the 
dust grains at radius r where the distribution of the maximum 
temperature T max (r) is steeper 1 (decreasing more quickly with 
radius) than considered for the clumpy case above. Therefore, 
the normalized surface density distribution, {rlR m ) a , has to be 
shallower, closer to the upper limit of a — 0. If there is a geo- 
metrical flaring, then we would revert to consider the directly il- 
luminated surface at each radius, leading to essentially the same 
conclusions obtained for the clumpy case. 

As for the consistency with the near-IR data, the large grain 
case produces slightly more compact distributions in the near- 
IR than the ISM case (because the temperature distribution is 
steeper), and is marginally consistent with the NGC4151 data. 
Thus the large grain case is slightly favored here. We note that 
the large grain case is consistent with the small innermost radii 
suggested by the near-IR rev erberations, though othe r contribut- 
ing factors are not ruled out ( IKishimoto et all [20071) . While we 
would need a more detailed modeling of the surface brightness 
distribution for the innermost region, we also need more obser- 
vational data in the near-IR to constrain such detailed models. 



(Bar vainii, Il987l) . For the large grain limit, y—0, while for the 
standard interstellar material (ISM) dust grains, y ~ 1.6. The 
sublimation temperature r su b is taken here as 1500K, but our 
conclusions below are not sensitive to this exact value of r su t,. 

The specific intensity I v from an illuminated dust cloud 
needs to be obtained with a proper radiative transfer calculation, 
but is approximately of the form 

I Y = J e- T B v (T(T)) dr (6) 

for an optically thick cloud, where T(r) is the temperature of 
the dust grains at an optical depth r from the surface of the 
cloud, and B V (T) is the Planck function. At a given observ- 
ing frequency, however, we can assume that the intensity dis- 
tribution over the torus radius, 7 y (r), approximately scales with 
B v (T m3X (r)). 

The second term in the expression for S v , (f/R m ) a , describes 
the distribution of the normalized surface density, or the filling 
factor per unit area, of the illuminated material as a power law 
form. An implicit assumption here is that the dust cloud distribu- 
tion in the innermost region is opaque so that S v (Ri n ) = / v (T SU b), 
and this is the maximum brightness so that a < 0. 

In this way, we can model the visibilities as a function of 
the parameter a for an assumed dust grain size distribution fixed 
by y. We consider the two cases above; namely, the large grain 
limit (y = 0) and the ISM case (y = 1.6). For each case, we can 
estimate an adequate density index a for roughly reproducing 
the visibilities observed in the 8-13 fim range. We show in Fig[2] 
such representative curves where a is taken to be 0.0 and -0.9 for 
the y = and 1.6 cases, respectively. We estimate the uncertain- 
ties in these indices to be ~0.3, based on the uncertainty in the 
visibility measurements and also in R\ n (~0.2 dex; see sec. 12. 11 ). 

All the curves shown in Fig[2]have been calculated with an 
outer boundary of 100 /?;„. The results do not change signifi- 



4. Conclusions 

The long-baseline IR interferometry data for AGN tori are still 
generally limited in uv coverage for each target. The data for var- 
ious objects probe different spatial scales. We argue that one way 
to uniformly study these various data is to view them as a func- 
tion of spatial frequency per inner torus boundary radius ^; n , or 
spatial scale in units of R m . In this way, using primarily the data 
for Type 1 AGNs, we have tried to construct a composite visibil- 
ity function over a wide range of spatial scales and investigated 
the radial structure of the tori. The data obtained so far suggest a 
common radial distribution of the face-on surface brightness in 
the mid-IR that is approximately of a power-low form with in- 
dex 2 and extends to ~100/?j n . Considering the temperature 

distribution of the dust grains, this corresponds to the surface 
density distribution of the directly illuminated material ranging 
approximately between r° and r~ l . We aimed here to derive di- 
rect constraints with a simple model. Further data are definitely 
needed over the wavelengths from the near-IR to mid-IR to test 
the composite visibility function and to further constrain models. 
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